AD-A071  856 


AERONAUTICAL  RESEARCH  ASSOCIATES  OF  PRINCETON  INC  NJ  F/6  20/4 

PRELIMINARY  STUOY  BY  MEANS  OF  SECOND-ORDER  CLOSURE  OF  THE  HEAT  —ETC(U) 
JAN  78  C D DONALDSON*  A K VARNA  N60921-77-C-0242 


n 


■ • ' *•  ••"**  - 

\T:&e  ■ 

- ■ 


PRELIMINARY  STUDY  BY  MEANS  OF 
SECOND-ORDER  CLOSURE  OF  THE  HEAT 
TRANSFER  TO  A POROUS  SPHERE/CONE 

EFFECTS  OF  SURFACE  ROUGHNESS, 
EDGE  TURBULENCE  AND  TRANSPIRATION 


NAUTICAL  RESEARCH  ASSOCIATES  OF  PRINCETON,  INC. 
i.  50  WASHINGTON  ROAD,  P.O.BOX  2229 
PRINCETON,  NEW  JERSEY  08540 


JANUARY  1978 

r ■y~:  ".i  Vr *■’ . ••  • 

ORT  FOR  PERIOD  30  SEPTEMBER  1977-  30  JANUARY  1978 


PREPARED  FOR 

NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK  LAB 
SPRING,  MARYLAND  20910 


SECURITY  CLktt'flCATION  or  TNI«  »A0l  i'S’HYI  Q«l»  t'm ered) 


REPORT  DOCUMENTATION  PAGE 


UEPC-Ki:  CCAJP1.ET1NC,  KOKM 


\l  OO VT  ACCfSJICN  MO 


TYY»e-qP  REJURT  l-VERIOD  COVERED 

Final  ifepSt*  ^ 

30  Sepfe*2»77-3l  Jan*4*78 


flrtf  bubrtrto). — 

PRELIMINARY  STUDY  BY  MEANS  OF  J5EC0ND-0RDER  Cl 
OF  THE  HEAT  TRANSFER  TO  A POROUS  SPHERE/CONE 

Effects  of  Surface  Roughness,  Edge  Turbulenci 

: *> 


A.R.A.P.  Report  No.  329 


I CONI  N ACT  0*  GRANT  KUMIC  R(-«J 


Ashok  K ./  Varma 


Coleman  duP 


Donaldson 


to.  program  element,  project,  task 


• PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Aeronautical  Research  Associates  of  Princeton 
50  Washington  Road,  P.O.  Box  2229 
Princeton,  NJ  08540 


APE  A A WORK  UNIT  NUMBERS 

61153N;  SF32-322 ; 

18635;  SEA  18451-SUM-32 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Surface  Weapons  Center,  WOL  / / I 

White  Oak,  Silver  Spring,  MD  20910 

Attn:  Code  WA-43  

U.  MONITORING  AGENCY  NAME  S AOORESSd/  dlllerent  Irom. Controlling  Qtlleo) 


IS.  SECURITY  CLASS.  (ol  title  r.porfj 


UNCLASSIFIED 


IS*.  DECL  ASSIFICATION  OOWNGRAOING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  ( ot  thle  Report) 


Approved  for  public  release.  Distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (ol  Mi*  akairacl  entered  In  Black  JO,  II  dlllerent  hem  Report) 


IS.  KEY  WORDS  (Continue  on  reeeree  elde  II  neeeeea ry  end  Identity  by  block  number) 

Turbulent  Boundary  Layers  Transition 

Blunt  Bodies  Roughness 

Conical  Bodies  Heat  Transfer 

Second-Order  Closure  Modeling 
of  Turbulence 

>0.  'Abstract  (Continue  on  ravarcc  elde  II  nacaaaary  end  Identify  by  block  number) 

nJA  compressible  second-order  closure  code  for  calculating  the  behavior  of 
boundary  layers  on  blunt  axisymmetric  bodies  at  zero  angle  of  attack  has  been 
The  code  is  capable  of  calculating  these  boundary  layers  both  with 


developed 

and  without  rotation  and  both  with  and  without  surface  mass  injection.  Pre- 
liminary calculations  have  been  made  with  the  code  to  compare  with  NSWC  measure 
iwHits  on 


sphere /cone  bod 
8 * IQ*;  and  vJJ.,6  * 10? 


M ■ 5 and  Reynolds  numbers  per  foot  of 
The  effects  of  surface  roughness  and  edge 
The  initial  theoretical  predictions 


turbulence  have  been  investigated 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THII 


(#>»•«  Dmf 


I 


UNCLASSIFIED 


HCU»‘T'  CLASSIFICATION  ©F  THIS  FAOtfHTiwi  C»i«  £nl»r»rf) 


20.  ABSTRACT  (concluded) 

rom  the  code  are  in  good  agreement  with  the  experimental  results.  The  results 
are  dependent  on  the  choices  for  the  surface  roughness  and  the  edge  turbulence 
parameters  and  it  appears  that  a more  extensive  parametric  search  could  lead 
to  even  better  agreement  with  the  experimental  data, 

' \ 


ICCUftlTv  CLASSIFICATION  OF  THIS  NAOeOFAm  Data 


1.  INTRODUCTION 


Over  the  past  few  months,  under  funding  from  the  Naval 
Surface  Weapons  Center  at  White  Oak,  A.R.A.P.  has  made  operational 
a compressible  second-order  closure  code  for  computing  the  be- 
havior of  the  boundary  layer  on  blunt  axisymmetric  bodies  at 
supersonic  Mach  numbers  when  the  angle  of  attack  is  zero.  The 
code  is  capable  of  calculating  such  boundary  layers  both  with  and 
without  rotation  and  both  with  and  without  continuous  surface  mass 
injection  (transpiration)  of  a gas  similar  to  that  in  which  the 
body  is  immersed.  In  the  process  of  checking  out  the  program, 
calculations  were  made  of  the  behavior  of  the  boundary  layer  and 
the  turbulent  heat  transfer  to  the  surface  of  a porous  sphere/cone 
without  rotation  but  both  with  and  without  transpiration.  The 
configuration  for  which  calculations  were  made  was  that  of  the 
sohere/cone  tested  by  Mr.  R.  H.  Feldhuhn  at  NSWC  at  Mach  number  5 
(Ref.  1).  The  Mach  number  and  Reynolds  numbers  for  which  calcu- 
lations were  made  were  those  reported  in  the  reference,  namely, 

M =5  and  ReM  = u R../N  * 6.33  x 10^  and  Re..  * 2.95  x 10^  . 
Here  R^  is  the  radius  of  the  SDherical  portion  of  the  body. 

For  these  conditions  parametric  studies  of  the  effects  of  surface 
roughness  and  boundary-layer-edge  turbulence  were  made  to  see  if 
the  code  could  predict  the  reported  results  for  reasonable  values 
of  these  two  parameters.  In  this  way  it  was  hoped  that  an 
evaluation  of  the  usefulness  and  capability  of  the  new  code  to 
predict  boundary  layer  phenomena  on  reentry  nose  tips  might  be  made 


2 . THE  CODE 


The  basic  equations  that  are  used  to  predict  compressible 
boundary  layer  behavior  on  blunt  bodies  of  revolution  are  given 
in  Ref.  2.  Also  given  there  are  the  details  of  the  modeling  that 
has  been  used  to  obtain  closure  of  the  original  set  of  equations. 
Values  for  the  coefficients  that  appear  in  this  compressible  model 
have  not  been  chosen  nor  are  they  adjusted  in  any  way  to  fit  the 
particular  set  of  data  that  is  considered  here.  Those  of  the 
coefficients  that  appear  also  in  incompressible  flow  modeling  are 
set  equal  to  the  values  that  have  given  the  best  fit  to  a large 
number  of  incompressible  flows,  many  of  which  are  not  boundary 
layer  flows.  The  rest  are  evaluated  by  analogy  with  their 
counterparts  in  incompressible  flow,  or,  where  no  counterpart 
exists,  they  are  set  to  zero.  More  than  a quarter  of  the  4l  co- 
efficients are  set  equal  to  zero  in  this  way.  In  the  sense  that 
the  same  parameters  are  used  for  all  flows,  the  full  model  is 
said  to  be  an  invariant  model  of  shear  flows  in  which  turbulence 
plays  a part.  In  the  case  of  no  rotation,  the  complete  model 
consists  of  a closed  set  of  13  partial  differential  equations  for 
the  following  quantities: 

u(s,n)  The  mean  velocity  parallel  to  the  local  surface  of 

the  body.  The  equations  specify  this  and  all  other 
dependent  variables  as  a function  of  distance  along 
the  body  s and  the  normal  distance  from  the 
surface  of  the  body  n . 


v(s,n) 

T(s,n) 

p(s,n) 

u,2(s,n) 


The  mean  velocity  normal  to  the  local  surface. 

The  mean  value  of  the  local  temperature. 

The  mean  value  of  the  local  static  pressure,  assumed 
given  at  the  outer  edge  of  the  boundary  layer. 

The  mean  value  of  the  square  of  the  velocity  fluctu- 
ations parallel  to  the  local  surface  of  the  body. 
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v'2(s,n) 

w,2(s,n) 


u' v' (s,n) 


T,2(s,n) 


p'u' 


p’ V1 


T'u’ 


T*v« 


The  mean  value  of  the  square  of  the  velocity  fluc- 
tuations normal  to  the  local  surface  of  the  body. 

The  mean  value  of  the  square  of  the  velocity  fluc- 
tuations perpendicular  to  u'  and  v'  . 

The  mean  value  of  the  correlation  of  the  velocities 
parallel  and  normal  to  the  local  surface. 

The  mean  value  of  the  square  of  the  temperature 
fluctuations. 

The  mean  value  of  the  correlation  of  the  fluctu- 
ations In  density  with  the  fluctuations  in  velocity 
parallel  to  the  local  surface. 

The  mean  value  of  the  correlation  of  the  fluctu- 
ations in  density  with  the  fluctuations  in  velocity 
normal  to  the  local  surface. 

The  mean  value  of  the  correlation  of  the  fluctu- 
ations in  temperature  with  the  fluctuations  in 
velocity  parallel  to  the  local  surface. 

The  mean  value  of  the  correlation  of  the  fluctu- 
ations in  temperature  with  the  fluctuations  in 
velocity  normal  to  the  local  surface. 


In  addition  to  the  above  quantities  we  have  the  equations  of 
state  which  relate  the  mean  density  p to  the  mean  temperature 
T when  the  mean  pressure  is  known.  Thus 


p = R ( pT  + p ' T ' ) 


(2.1) 


(2.2) 


In  the  present  model,  as  indicated  above,  p’T*  is  modeled  as 
_(p/T)T'2  . 


The  code  permits  one  to  solve  the  set  of  coupled  partial 
differential  equations  for  the  above  variables  as  functions  of 
s and  n given  the  following  information: 

1.  The  shape  of  the  blunt  body  and  its  inviscid  pressure 
distribution . 

2.  The  initial  values  of  all  13  dependent  variables  as  a 
function  of  n at  some  station  which  may  be  close  to  the  stagna- 
tion point  but  such  that  s/6  >>  1 where  6 is  the  initial 
boundary  layer  thickness. 

3.  Conditions  at  the  edge  of  the  boundary  layer  and  the 
temperature  and  transpiration  wall  conditions,  all  as  functions 
of  s . 

For  the  runs  discussed  here,  it  is  assumed  that  the  flow  at 
the  initial  station  is  close  to  the  solution  for  a completely 
laminar  stagnation  point.  It  is  possible  to  remove  this  restric- 
tion, but  a program  is  not  yet  available  for  solving  for  a stag- 
nation point  boundary  layer  in  the  presence  of  free-stream 
turbulence.  The  development  of  such  an  initializing  code  would 
seem  a logical  complement  to  the  existing  program  if,  indeed, 
this  program  can  be  shown  to  give  insights  into  the  problem  of 
transition  on  blunt  bodies. 

Before  going  on  to  discuss  some  preliminary  calculations  of 
blunt  body  heat  transfer  that  were  made  to  test  the  code,  mention 
should  be  made  of  the  way  in  which  the  turbulence  scales  were  set 
in  these  calculations.  When  using  A. R. A. P.’s  second-order-closure 
codes  to  predict  turbulent  phenomena,  it  is  always  necessary  to 
specify  or  calculate  the  local  scale  of  the  turbulent  fluctuations 
A(s,n)  . A differential  equation  for  A may  be  written,  but  the 
coefficients  of  the  terms  in  this  equation  depend  on  the  anisot- 
ropy of  the  turbulence  that  is  produced  and  are  not,  therefore, 
invariant.  In  view  of  this,  we  have  chosen,  for  the  present,  to 
specify  the  scale  A(s,n)  as  a function  of  certain  parameters 
of  the  mean  flow  in  a way  that  has  been  shown  to  give  good 
results  for  turbulent  boundary  layers.  For  this  reason,  A(s,n) 
is  taken  to  be  the  minimum  of  the  two  quantities  A^n  and 


A where 
out 

* Atn'Ao+0'65n  (2.3) 

Aout  ’ 017S.99<6>  <2A> 

c 

The  parameter  A is  used  to  represent  the  effects  of 
surface  roughness.  It  has  been  found  empirically  that  AQ  is 
approximately  proportional  to  the  magnitude  of  the  surface 
roughness.  Indeed,  it  has  been  found  that  a useful  equivalence 
between  surface  roughness  of  height  k and  the  parameter  Aq 
is 

6.5Ao  < k < 13Aq  (2.5) 

£ 


The  first  use  made  of  the  A.R.A.P.  blunt  body  boundary  layer 
code  was,  as  mentioned  in  the  introduction,  the  prediction  of  the 
behavior  of  the  boundary  layer  on  the  porous  sphere/cone  tran- 
spiration cooled  model  tested  at  M = 5 at  NSWC/White  Oak. 


Initial  conditions.  For  the  calculations  reported  here  the 
following  initial  conditions  were  specified  at  s/RN  = Sj/RN  = 0,1  * 
(The  corresponding  values  of  s/6  were  in  the  range  30  < s/6  < 70  .) 
The  boundary  layer  was  assumed  to  be  of  a form  close  to  the  appro- 
priate laminar  stagnation  point  solution.  It  was  also  assumed  that 
a small  velocity  disturbance  existed  in  this  initial  boundary.  In 
particular,  for  the  calculations  reported  here,  the  velocity  dis- 
turbance was  assumed  isotropic  in  the  sense  that 


(u'u'Jj  * (v'v')1  = (w'w')i  = q^/3 


(3-1) 


The  distribution  of  this  initial  turbulent  energy  was  taken  to  be: 

2 

q1(s1,n)  = 0 for  0 < n < 0.126  ^(s^)  j from  n = 0.126  gg(s^) 
to  n = 0.26  gg(s1)  , q^(s1,n)  was  linear  in  n and  reached  its 
maximum  value  q?  = 3 x 1 0~5  u2  at  n = 0.26  QQ(s.  ) ; from 
n = 0.26  to  n = 0.35  99^1^  > q (s^n)  decreased  linearly 

with  n and  reached  zero  at  n * 0. 35  9 9 ( s ^ ) ; for 
n > 0.36  ^(s^)  » q^  was  assumed  equal  to  zero.  All  other  turbu- 
lent correlations  were  assumed  zero  at  si  . 

Calculations  for  zero  roughness.  In  Figs.  3-1  through  3.^ 
we  show  calculations  of  the  heat  transfer  coefficient  for  the 
test  sphere/cone  of  Ref.  1 at  Reynolds  numbers,  Re^  = u^Rj^/N^  , 
of  6,3  x 10^  and  2.95  x 10^  both  with  and  without  surface 
transpiration.  The  cases  calculated  correspond  to  cases  1,  2,  22, 
and  81  of  Ref.  1.  The  calculations  shown  were  made  assuming  zero 
surface  roughness  (Aq  = 0)  and  zero  boundary-layer-edge  turbulence. 
For  both  high  and  low  Reynolds  numbers  where  there  is  no  transpi- 
ration, the  calculations  indicate  that  the  boundary  layer  will  not 
undergo  early  transition  to  turbulent  flow  for  the  small  initial 
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Ref.  1,  Case  no.  22 
With  transpiration 

A.R.A.P.  predictions 
A0S  0 ks  0 
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O Experimental  data 


Figure  3.4.  Calculated  heat  transfer  coefficient  for  conditions 
of  Case  No.  22  (Ref.  1)  for  zero  surface  roughness  and  zero  edge 
turbulence  . = A 'X'X  x Ifi-)  wf  hVi  hrencnfrofinn 


6.33  * 1CP  with  transpiration. 


spot  of  turbulence  we  have  assumed  (Figs.  3.1  and  3.2).  However, 
when  the  destabilizing  effect  of  transpiration  is  included  in  the 
calculations,  we  find  that  the  higher  Reynolds  number  case  ex- 
hibits early  transition  while  the  low  Reynolds  number  case  does 
not  (Figs.  3*3  and  3-4).  We  note  that  the  calculated  heat  trans- 
fers for  the  high  Reynolds  number  case  when  there  is  zero  rough- 
ness and  zero  free-stream  turbulence  are  below  those  measured  on 
the  test  sphere/cone.  It  is  also  clear  from  a detailed  study  of 
the  computer  output  that  the  whole  behavior  of  the  flow  in  the 
stagnation  point  region  is  extremely  sensitive  to  the  turbulence 
that  is  introduced  into  the  boundary  layer  in  this  region. 

Indeed,  the  occurrence  of  early  transition  will  depend  on  a com- 
plex interaction  of  the  Reynolds  number,  the  free-stream  turbulence 
level,  the  body  roughness,  the  body  temperature,  the  Mach  number, 
and  so  forth. 

Calculations  with  surface  roughness.  In  Figs.  3.5  and  3.6 
we  show  the  effect  of  surface  roughness  for  the  case  of  high 
Reynolds  number.  For  the  case  when  there  is  no  transpiration 
(Case  1),  we  show  in  Fig.  3.5  the  effect  of  roughness  through 
calculations  having  uniform  roughness  over  the  entire  sphere/cone 
of  the  following  levels: 


A = 
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k = 0 


= 3 * 

10-6 

inches 

20 

< 

k 

< 

40  microinches 

= 3 x 

10"5 

inches 

200 

< 

k 

< 

400  microinches 

= 3 x 

«=T 

1 

O 

rH 

inches 

.002 

< 

k 

< 

.004  inches 

It  Is  clear  from  Fig.  3*5  that  a surface  roughness  in  the 

range  from  20  to  40  microinches  will  not,  according  to  these 

particular  calculations,  cause  early  transition  at 

Re,.  = 2.95  x 10^  when  there  is  no  steady  source  of  edge  turbu- 
N 

lence.  However,  if  the  roughness  is  increased  to  the  order  of 
from  200  to  400  microinches,  early  transition  will  occur.  The 
level  of  heat  transfer  In  the  nose  region  is,  by  calculations, 
somewhat  less  than  is  observed  experimentally  for  this  level  of 
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Figure  3.6.  Effect  of  two  surface  roughnesses  on  heat  transfer 
with  edge  turbulence  equal  zero  for  Case  No.  81  (Ref.  1). 

ReN  - 2.95  x 106  with  transpiration. 


roughness.  For  roughness  in  the  range  from  .002  to  . 0044  inches, 
early  transition  occurs  and  the  levels  of  heat  transfer  in  the 
nose  region  are  more  nearly  equal  to  the  experimental  values.  For 
this  larger  roughness  the  calculated  heat  transfer  to  the  conical 
region  of  the  body  is  considerably  higher  than  the  experimental 
data.  Although  not  shown  in  Fig.  3.5,  calculations  carried  out 
for  roughness  of  the  order  of  .001  to  .002  inches  fell  almost  on 
top  of  the  results  for  roughness  of  the  order  of  .002  to  . OOh 
inches  on  the  spherical  portion  of  the  body,  and  between  the 
curves  for  200  < k < *400  microinches  and  .002  < k < .00*4  inches 
on  the  conical  section  of  the  body. 

For  the  case  when  transpiration  is  considered  at  high 
Reynolds  number  (Case  81)  we  show  calculations  in  Fig.  3.6  for 

A - 0 k - 0 

o 

-44 

A ■ 3 x 10  inches  .002  < k < . 0044  inches 

o 

Here  again  we  see  that  roughness  Increases  the  turbulent  heat 
transfer  in  the  stagnation  point  region  to  values  close  to  the 
experimental  data  on  the  spherical  portion  of  the  body  as  well  as 
on  the  cone.  As  before,  for  levels  of  roughness  that  give  heat 
transfer  in  agreement  with  experimental  results  on  the  spherical 
portion  of  the  body,  the  heat  transfer  to  the  conical  portion  of 
the  body  is  overpredicted.  It  appears  that  we  are  dealing  with  a 
problem  where  both  the  roughness  level  and  boundary-layer  edge 
turbulence  level  are  needed  to  achieve  a match  to  experimental 
data  in  both  regions. 

Calculations  with  roughness  and  free-stream  turbulence.  The 
final  calculations  we  shall  show  to  round  out  this  preliminary 
study  of  the  A.R.A.P.  code  are  calculations  where  we  have  fixed 
the  surface  roughness  at  a value  guessed  at  from  the  data  Just 
presented  and  varied  the  level  of  external  turbulence  impinging  on 
the  edge  of  the  boundary  layer.  The  same  initial  conditions  as 
those  used  in  the  previous  runs  were  used.  This  is  not  exactly 
correct,  for  if  there  were  external  turbulence,  the  turbulence 
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distribution  at  s/R^  *0.1  would  not  be  a spike  in  the  interior 
of  the  boundary  layer  with  a uniform  outer  turbulence.  However, 
in  making  these  first  calculations  some  parameters  of  the  problem 
had  to  be  held  fixed. 

The  level  of  roughness  chosen  was  .001  < k < .002  inches. 
This  roughness  was  placed  only  where  there  were  porous  surfaces 
on  the  model,  that  is  for  s/R^  < 3.24.  The  roughness  was  taken 
to  be  zero  outside  of  this  region.  Two  levels  of  boundary-layer- 
edge  turbulence  were  investigated.  They  were  q /xr  * 3 * 10” 

— 7 0 00 

and  3 x 10”  . This  turbulence  was  assumed  to  be  isotropic  in 

the  sense  that  u'2/u2  » v'2/u2  * w'2/u2  . The  results  are  shown 

0 OO  0 OO  0 00 

in  Figs.  3.7  through  3.10. 

For  the  high  Reynolds  number  cases  (Figs.  3.7  and  3-8),  it  is 

2 -2 

not  worth  plotting  the  curves  for  the  separate  values  of  q /u  , 

0 OO 

because  the  curves  are  not  markedly  different.  In  these  cases  the 
curves  differ  only  in  a small  region  close  to  the  stagnation  point 
and,  without  a proper  initializing  code,  the  differences  in  the 
results  are  probably  not  too  meaningful.  However,  for  the  low 
Reynolds  number  cases  (Figs.  3-9  and  3-10)  the  effect  of  free- 
stream  turbulence  is  pronounced.  Transition  is  occurring  on  the 
conical  portion  of  the  test  body  for  these  cases  and  the  computed 
results  are  not  in  bad  agreement  with  the  experimental  data.  The 
rate  at  which  transition  is  occurring  does  appear  to  be  too  fast. 
However,  this  may  be  due  to  the  distribution  of  edge  turbulence 
that  has  been  assumed. 

Before  going  on  to  make  some  general  conclusions  and  recom- 
mendations for  what  might  be  future  work,  the  effect  of  cutting 
off  the  roughness  for  the  high  Reynolds  number  cases  at  the  end 
of  the  porous  wall  region,  s/R^  " 3*2^  should  be  noted.  The 
rapid  drop  in  heat  transfer  at  this  point  and  the  fact  that  compu- 
tations give  too  high  a heat  transfer  on  the  cone  ahead  of 
s/Rn  ■ 3*2^  , would  suggest  that  the  roughness  we  have  chosen  is 
too  high  and  the  actual  porous  material  roughness  must  be  less 
than  .001  to  .002  inches.  This  deduction  is  strengthened  when  we 
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Figure  3.7.  Calculated  heat  transfer  coefficient  for  a particular 
choice  of  surface  roughness  and  edge  turbulence  for  Case  No . 1 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  given  a very  brief  description  of  the  results  ob- 
tained in  an  initial  study  of  the  effects  of  surface  roughness 
and  boundary-layer-edge  turbulence  on  the  heat  transfer  to  a 
blunt  body  with  and  without  surface  transpiration  using  the 
A.R.A.P.  second-order-closure  code.  It  seems  that  this  code 
exhibits  considerable  promise  for  the  calculation  of  boundary 
layer  behavior  on  high-speed  blunt  bodies.  Indeed,  it  would 
appear  that,  given  a more  complete  parameter  search,  a choice 
of  surface  roughness  and  edge  turbulence  could  be  found  that 
would  give  considerably  better  agreement  between  calculations 
and  the  experimental  data  we  have  so  far  examined.  The  values 
of  roughness  and  turbulence  that  would  do  this  appear  to  be  very 
realistic . 

It  is  recommended  that  a more  extensive  parametric  study  of 
the  effects  of  surface  roughness  and  edge  turbulence  than  that 
carried  out  in  demonstrating  this  code  be  carried  forward  using 
the  code  in  Its  present  form. 

It  is  clear  that  an  initializing  code  which  will  allow  proper 
computation  of  stagnation  point  flows  when  the  impinging  flow 
contains  turbulence  fluctuations  is  desirable. 

It  is  recommended  that  such  an  initializing  code  be 
developed  forthwith. 
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